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THE  BOUN0ED  ^-POINT  CIASSES  IN  ROBUST  HYPOTHESIS  TESTING  AND  FILTERING* 
SALEEM  A.  KASSAM 

Department  of  Systems  Engineering 
Moore  School  of  Electrical  Engineering 
University  of  Pennsylvania 
Philadelphia,  PA.  19104 

ABSTRACT 

The  bounded  p-point  classes  of  density  functions  are  defined  through 
upper  and  lower  bounds  on  the  power  or  probabilities  on  different  sets  in 
some  partition  of  the  frequency  or  sample  space.  The  robust  Wiener  filter 
is  derived  for  signal  and  noise  spectra  lying  in  such  classes,  and  is 
shown  to  be  piecewise  constant  in  its  frequency  response.  Similar  results 
are  indicated  for  hypothesis  testing  involving  such  classes  of  probability 
densities.  The  robust  matched  filter  is  also  obtained  for  the  bounded  p- 
point  class  of  noise  spectra,  and  is  again  shown  to  have  a  piecewise- 
constant  frequency  response. 

I.  INTRODUCTION 


In  recent  years  a  considerable  number  of  results  has  been  obtained  on 
mlnimax  robust  filters  operating  under  conditions  of  inexact  a  priori 
knowledge  of  signal  and  noise  characteristics.  Specifically,  filters  for 
estimation,  smoothing,  prediction,  and  interpolation  of  random  signals 
have  been  considered  [e.g.,  1-3],  in  addition  to  matched  filters  for 
signal-to-noise  ratio  optimization  (e.g.,  4],  and  other  filters  for  special 
applications,  e.g.,  time-delay  estimation.  In  all  these  cases  a  priori 
uncertainties  on  signal  and  noise  characteristics  are  modeled  by  classes 
of  possible  second-order  characteristics,  i.*".,  spectral  densities  or  co- 
variance  functions.  Many  of  these  results  were  motivated  by  earlier  work 
on  robust  hypothesis  testing  [5],  in  which  minimax  robust  tests  for  binary 
hypotheses  were  obtained  for  some  specific  classes  of  probability  density 
functions.  In  particular,  minimax  robust  solutions  for  robust  Wiener  fil¬ 
tering  under  spectral  uncertainty  classes  and  for  robust  hypothesis  testing 
under  corresponding  probability  density  function  classes  are  very  closely 
related. 

Although  some  rather  general  results  are  available  on  robust  hypothesis 
testing  and  robust  Wiener  filtering  under  general  uncertainty  classes,  sol¬ 
utions  for  useful  specific  classes  are  of  particular  interest  in  applications. 
The  spectral  and  probability  density  models  for  which  specific  solutions 
have  been  obtained  include  the  following:  1)  the  e-contaminated  model  [5]; 

2)  the  "band-model"  for  bounded  densities  [6];  3)  the  bounded  total-var¬ 
iation  model  [5];  and  4)  the  specified  p-point  model  [7]. 


Of  these,  perhaps  of  primary  Importance  in  applications  are  classes 
generated  by  the  "band-models"  and  the  specified  p-point  models.  In  the 
band-classes  density  functions  are  assumed  to  lie  between  upper  and  lower 
bounding  functions,  and  in  the  specified  p-point  classes  density  functions 
are  assumed  to  have  specified  proportions  of  their  total  area  in  specified 
regions  in  the  sample-space  or  frequency  domain.  A  natural  and  useful 
synthesis  of  these  two  types  of  classes  yields  the  bounded  p-point  class 
of  densities,  where  density  functions  are  assumed  to  have  proportions  of 
their  total  area  in  specified  regions  lying  between  specified  upper  and 
lower  bounds.  It  is  easily  seen  that  the  bounded  p-point  classes  include  as 
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a  limiting  case  the  band-classes,  and  as  a  special  case  the  specified  p- 
point  classes. 

We  obtain  the  structure  of  robust  Wiener  filters  for  signal  and  noise 
spectra  in  bovinded  p-point  classes  defined  by  a  common  partition  of  the 
frequency  space,  by  using  general  results  on  robust  hypothesis  testing  for 
bounded  classes  [6] .  The  problem  is  first  converted  to  one  involving  dis¬ 
crete  spectra,  and  then  results  on  robust  hypothesis  testing  for  discrete 
probability  distributions  are  used.  The  robust  Wiener  filter  has  a  piece- 
wise  constant  frequency  response.  In  a  similar  way  it  can  be  shown  that 
robust  hypothesis  testing  for  bounded  p-point  classes  of  density  functions 
requires  the  use  of  a  quantizer  operating  on  input  data.  For  Wiener  fil¬ 
tering  we  also  obtain  a  result  for  the  worst-case  performance  of  any  filter 
when  signal  and  noise  spectra  lie  in  bounded  p-point  classes  or  the  bounded 
classes  of  the  "band"  model.  The  above  results  are  obtained  for  signal  and 
noise  classes  defined  on  the  same  partition  of  the  sample  or  frequency  space. 
We  also  show  that  a  solution  can  be  obtained  when  one  of  the  densities  is  in 
the  bounded  class  generated  by  the  band-model. 

For  matched  filtering  under  noise  spectral  density  uncertainty,  we 
also  consider  the  bounded  p-polnt  class  for  the  noise  spectra.  For  known 
signal  the  solution  is  simple.  When  the  signal  is  also  uncertain,  lying  in 
the  previously  used  bounded  L2-deviation  class  [4],  we  can  also  obtain  a 
general  solution.  The  interesting  result  here  is  that  the  robust  matched 
filter  frequency  response  is  again  a  plecewise-constant  function. 

II.  ROBUST  HYPOTHESIS  TESTING  AND  WIENER  FILTERING 

We  start  by  seeking  a  solution  for  the  robust  Wiener  filter  when  signal 
and  noise  power  spectral  densities  (PSD's)  are  members  of  bounded  p-point 
classes.  This  will  involve  the  robust  binary  hypothesis  test  for  discrete 
"band"  classes  of  probability  densities.  We  define  the  bounded  p-polnt  class 
S  of  allowable  signal  PSD's  S  by 

s  -  {S:  /.  S(f)df  =  o;:  ,  j  =  1.2,...,m,  I  oT  =  c^}  (la) 
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Here  ™  is  a  given  partition  of  the  frequency  space  (with  f  t  A. 

1  J  2  J 

-f  e  A  ),  and  is  the  known  total  signal  power,  with  the  fractional  powers 
2  ^ 

Oj  satisfying 

^  ■  ‘4j  ’  j  "  l,2,...m  (lb) 

2  2 

where  the  bounds  ^  ,  0^^  are  known.  In  a  similar  way  we  define  the  class 
W  of  allowable  noise  PSD's  N  by 


M  =  {N:  N(f)df  =  V,,  j  =  l,2,...m,  ^  v  =  v  > 

Aj  J  J-l  J 
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Note  that  we  are  assuming  here  a  common  partition  both  signal 

and  noise  classes. 


We  now  consider  a  reduced  discrete  spectral  distribution  problem.  Let 


all  signal,  noise  power  in  each  subset  A.  be  concentrated  at  a  single  fre¬ 
quency  pair  ±f.  e  A  .  We  will  denote  by^  5  ,  W  the  resulting  classes  of 
one-sided  power  spectral  concentrations  derIved“from  ,  W,  respectively. 

Let  H  be  the  frequency  response  of  any  filter  used  for  estimating  tlie  sig¬ 
nal  component  with  power  spectrum  =  [S  (1),  5^(2) , . . . S^Cm) ]  in  S 
observed  in  additive  uncorrelated  noise  with  power  spectrum  N  ^  W  .  Let 
Hj(j)  =  H(f^),  j  =  1,  2,...,m.  Then  the  mean-squared-error  Onse) 
^d^^d^d’^d^  in  estimating  the  signal  for  this  discrete  PSD's  case  is 

m  „ 

^d^^d*^d‘"d^  +  N^(j)iH^(j)i-  (3) 

For  given  S^,  this  is  minimized  for  H^(j)  =  S  ,(j)/[S  (j)+N  (j) ]  and  the 
corresponding  minimum  mse  is 

*  ^ 

^  S^(j)N^(j)/[S^(j)+N^(j)]  (4) 

j=l 

.1  I , 

To  obtain  the  least  favorable  pair  (S^.N’"),  which  is  defined  to  be  the 
pair  maximizing  (4),  we  examine  the  corresponding  hypothesis  testing  pro¬ 
blem  for  discrete  densities.  Consider  a  random  variable  with  a  finite  set 
of  possible  realizations  {l,2,...,m},  with  a  probability  mass  function  given 
by  either  We  assume  that  Pq  and  P  are  in  the 

bounded  classes  defined  by  ^ 


^0  “  ^^0  '  ^  ’  j  =  l,...,m}  (5a) 

Pf  »  {P^  ;  P^,^(j)  ^  P^(j)  i  P^^.(j)  ,  j  =  l,2...m}  (5b) 

»  I 

Now  the  pair  of  densities  (Pq,  P^)  in  these  classes  which  is  least-favorable 
(risk)  in  testing  ?f]  vs.  [6]  can  be  found  by  applying  the  general  result 
in  [6]  for  "band"  classes  of  bounded  probability  density  functions.  The 
results  in  [6]  are  applicable  for  both  absolutely  continuous  and  discrete 
probability  distributions.  Note  that  the  definition  of  the  pair  (P^.Pp 
which  is  lea§t-f avorable  (risk)  implies  that  a  corresponding  optimum  test 
for  pi  vs.  P^'  is  robust  for  testing  P«  vs.  P^  Op  r);ie  other  hand,  the  pre- 
vious'^def inirion  of  the  least-favorable  pair  does  not  automat icallv 

imply  robustness  of  the  corresponding  optimal  filter,  although  tor  our  case 
this  will  be  shown  to  be  true.  In  order  to  solve  for  the  least-favorable 
pair  (S^,N^')  from  the  results  in  [6]  applied  to  discrete  probabilitv  dis¬ 
tributions,  we  use  the  following  lemma  which  is  a  direct  counterpart  of 

Lemma  1  in  [2]  where  it  is  stated  for  the  absolutely  continuous  case: 

Lemma  1  Let  Pq  and  P.  be  classes  of  discrete  probability  mass  functions 
P  =  such  that  each  member  of  Pq^P^  has  positive  values  on  the 

common  set  {l,2,...m}.  If  Pq  and  P^  are  least-favorable  (risk)  for  Pg 
then 


(b) 


The  proof  of  this  lemma  follows  the  same  ideas  as  those  used  in  the 
proof  of  lemma  1  in  [2].  Although  this  lemma  can  be  extended  to  apply  to 
the  slightly  more  general  case,  we  will  instead  restrict  our  original  classes 
of  spectral  densities  by  requiring  that  the  lower  bounds  .  and  be 
positive  (except  in  the  case  where  both  and  v^.  are  zer^).  Now^we  trans¬ 
form  the  original  classes  S,  and  into  classes  Pg  and  P  of  probability  mass 
functions  P  and  P  respectively,  by  scaling  each  member  of  Sj  by  l/o"^  and 
scaling  eac§  member  of  N,  by  1/v^.  Thus,  the  P  c  P  are  defined  by 
P  (j)  =  S^(j)/o^-,  j=l,2,...,m.  Let  (P^  P^-)  be  Ihe  llast-favorable  (risk) 

pair  in  P  x  P  .  Now  e§  of  (4)  can  written  as 
s  n  u 


e,(P  ,P  )  =  y  - Lli) - 

dsn  T  2 


PnO) 


where  L(j)  =  P  (j)/P  (j).  The  function  x/(o^x+v")  is  a  concave  function  of 

X,  and  lemma  l^eives'^the  result  that  SV  =  a‘P^  and  N^  =  v‘^P^  form  the  least 

•-*  Q  s  o  n 

favorable  pair  maximizing  (4)  over  the  classes  S,  and  W,. 

d  d 

Once  the  least  favorable  pair  (S^,N^)  has  been  characterized,  it  re¬ 
mains  to  be  shown  that  the  optimum  filter  for  this  signal,  noise  spectra 
is  the  robust  filter  minimizing  the  maximum  mse  over  the  classes  S,  and 
This  result  is  also  a  direct  counterpart  of  Theorem  1  in  [2]  for  the  con¬ 
tinuous  case,  and  may  be  stated  as  follows: 

Lemma  2  Let  5^,  W,  be  convex  classes  of  power  spectral  concentrations  on  a 
finite  number  of  frequencies.  Then  (S^,N^)  c  ,  x  Wj  is  least-favorable 
for  Wiener  filtering  iff  (S^,N^)  and  H  ,  =  SV*^Sj+Np  form  a  saddlepoint 

for  the  mse  functional  e^(S^,N^;H^) ,  ’ 

Proof  It  is  clear  that  if  a  pair  (S^,N^)  and  its  optimal  filter  _  form 
a  saddlepoint  for  mse,  then  the  pair  is  least-favorable.  To  prove’^e  con¬ 
verse  we  note,  as  in  the  proof  in  [2],  that  the  functional  is 

concave  in  S^and  in  N  ,  and  that  S^,  are  assumed  to  be  convex.  From 

this  the  proof  proceeds  in  a  manner  similar  to  that  in  [2],  except  that  now 
we  use  the  mse  expressions  (3)  and  (4)  in  place  of  the  error  expressions  for 
the  continuous  case. 

We  thus  find  that  the  saddlepoint  robust  filter  for  the  bounded  classes 
Nj  of  discrete  spectra  can  be  found  as  the  optimal  filter  for  the  least- 
favorable  pair  in  S.  x  W^.  This  least  favorable  pair  is  obtainable  as  a 
scaled  version  of  the  least-favorable  (risk)  pair  for  testing  the  hypo¬ 
thesis  P  vs.  P_  where  P  and  P  are  derived  by  scaling  members  of  S,  and  hi, 
s  s  ^  ^ 

respectively.  We  now  finally  obtain  the  saddlepoint  robust  filter  for 

the  original  bounded  p-point  spectral  classes  S  and  M. 

Lemma  3  If  H,  is  the  saddlepoint  robust  filter  for  the  classes  Sd,N^  of 
discrete  spectia  derived  from  the  bounded  p-point  classes  S,W  defined  by 
(1)  and  (2),  then  the  filter  defined  by 

Hj(£)  -  „a).  f  e  Aj  (8) 

is  the  saddlepoint  robust  filter  for  spectral  classes  S  and  hi 


Proof  We  have,  for  S  e  5  and  N  e  hi. 


e(S,N;Hj^)  -  /S(f)  Il-Hj^(f)  +  N(f)lH  (f)i^df 


m 

1  lil-H 

j=l 


d,R 


(j)i 


-sjo) 


=  cCs^.nShj^) 

*  a  z 

=  e  (s\n^). 


where  S  ,  N  is  any  PSD's  In  S,N  respectively  for  which  S^(f)/N^(t) 

=  Sj(j)/N^(j)  when  f  £  Aj .  Note  that  such  a  pair  is  least-favorable. 

What  we  have  obtained  is  an  explicit  solution  for  the  piecewise-constant 
robust  Wiener  filter  for  bounded  p-point  classes  of  signal  and  noise 
spectra.  In  a  similar  way,  we  can  obtain  a  corresponding  result  for  ro¬ 
bust  hypothesis  testing  for  bounded  p-point  classes  of  probability  density 
functions . 


Returning  to  the  Wiener  filtering  problem,  let  S  be  a  bounded  p..point 
class  of  signal  spectra,  defined  with  resp^t  to  some  partition  Let 

aJ  be  any  class  of  noise  spectra,  and  let  W  be  the  class^^of  bourid^d "Ap¬ 
point  spectra  generated  by  R,  for  the  same  partition  {A.}j=l.  We_  then  con¬ 
sider  the  least-favorable  pairs  (S^,N^)  for  these  classes  ^j^and  R  of  signal 
and  noise  spectra.  If  a  least-favorable  pair  exists  with  N  £  R  ,  then 
that  pair  is  also  least  favorable  for  the  original  passes  5  and  R  .  We  can 
apply  this  observation  to  the  case  where  the  class  W  is  a  "band"  rlass_^f 
bounded  spectra  as  in  [1].  From  the  proof  of  lemma  3  we  note  that  S^,N 
are  required  to  be  constant  multiples  of  each  other,  their  exact  shapes  not 
being  con£|rained.  Therefore  it_is  always  p£ssible  to  find  a  least-favorable 
pair  (S^,N  )  for  classes  S  and  N  such  that  is  in  the  band  cl^ass  R.  This 
implies  that  the  more  precise  bounds  on  noise  spectra  defining  W  represents 
extra  information,  beyond  bounds  on  fractional  powers,  which  is  not  useful 
in  obtaining  a  better  robust  filter.  Note  that  in  the  above,  the  roles  of 
S  and  R  can  be  interchanged  and  similar  statements  can  be  made.  For  hypo¬ 
thesis  testing,  again  similar  conclusions  can  be  drawn. 


Before  considering  matched  filters  in  the  next  section,  we  establish 
a  result  which  allows  computation  of  the  worst-case  performance  of  any  fil¬ 
ter  for  signal  and  noise  spectra  defined  either  by  "band"  models  or  the 
bounded  p-polnt  classes. 


In  general,  for  given  signal  and  noise  classes  it  is  of  interest  to 
have  the  worst-case  performance  of  any  given  filter  H.  Consider  the  power- 
constrained  band  classes  of  spectra  for  Wiener  filtering.  Note  that  the 
mse  obtained  with  any  filter  frequency  response  H  is  given  by 

e(S,N;H)  -  /N (f ) lH(f ) j ^df  +  /S (f ) | 1-H(f ) j ^df 


Consider  any  Integral  of  the  form  1  ■  jX(f)p(f)df  where  P(f)  is  given  and 
X(f)  lies  in  a  power  constrained  "band"  class  of  power  spectra.  Then  the 
worst  (maximum)  value  of  I  is  obtained  with  X  in  the  "band"  class  defined  by 


P(f)  <  k 
P(f)  >  k 
P(f)  =  k 


(10) 


fx.  (f). 

Arbitrary, 


where  k  is  a  constant  chosen  to  satisfy  the  power  constraint.  The  proof 
follows  easily  by  considering  /[X(f)-X  (f)]P(f)df  for  any  other  X.  This 
result  can  be  applied  to  each  of  the  two  terms  in  the  above  error  expression 
to  obtain  the  worst  case  signal  and  noise  spectra  in  the  "band"  class  for 
any  given  filter.  A  direct  extension  gives  the  worst-case  spectra  for  the 
bounded  p-point  classes. 


III.  BOUNDED  p-POINT  NOISE  CLASS  IN  ROBUST  MATCHED  FILTERING 


For  noise  characteristics  described  in  the  frequency  domain,  the  fre¬ 
quency  response  of  the  robust  matched  filter  was  derived  explicitly  in  [4] 
for  noise  spectral  uncertainty  modeled  by  the  "band"  classes  of  allowable 
spectra.  The  signals  were  required  to  be  within  a  given  L^  distance  of  a 
nominal  signal.  Here  we  will  obtain  results  for  the  case  where  the  noise 
class  is  the  bounded  p-point  class  described  by  (2).  If  S  is  the  Fourier 
transform  of  the  real  deterministic  signal  and  S  is  the  nominal  charac¬ 
teristic,  the  signal  class  S  is  now  defined  by 

S  =  {S:  /^ls(f)  -  S  (f)l^df  ^  6}  (11) 

Ji  o' 

Here  Q  is  (-«>,“>)  for  the  continuous  case  or  [-O' 5,  O' 5)  for  the  discrete 
case.  The  signal-to-noise  ratio  obtained  when  using  a  filter  frequency  re¬ 
sponse  H,  when  S  and  N  are  the  given  characteristics,  is  given  by 


SNR(S,N;H)  =  1 1 
^*N(f)|H(f)|'^df 


(12) 


where  integrals  are  over  Q.  The  following  theorem  gives  our  general  result 
defining  the  least-favorable  pair  (S  ,N  ) ,  which  together  with  the  corres¬ 
ponding  optimum  filter  H^  forms  a  saddlepoint  for  the  SNR  functional. 


Theorem  1  For  the  bounded  p-ooint  noise  class  of  (2)  and  the  bounded  L^- 

deviatlon  signal  class  of  (ll),the  pair  (Sj^.N^^)  defined  below  is  least- 

favorable  for  matched  filtering  and  the  optimum  filter  H^^  for  (Sj^,Nj^)  is  a 

saddlepoint  robust  (matched)  filter,  provided  non-negative  constants  k  and 

c  exist  satisfying  (17)  and  (18).  The  pair  (S  ,N  )  is  defined  by 

R  K 


=Rj">  - 


|(|S„(f)|  -  |S^(f)|>k.c  (13) 


,  Otherwise 


and 


NRj(f)  =  |Sj^j(f;[/k^. 


(14) 


for  j”l,2 . .  where  Sj^ . ,  N^^.  are  the  restrictions  of  Sj^,  to  the  sets 

A,,  respectively.  Here  R-]  are^defined  as  functions  of  N^^,  k  and  c 
through  the  equations  ^ 


4  ^ 


(15) 


and 


the  set  R.»  j  =  being  defined  as  that  subset  of  A  on  which 

!  ^  ^ 

S  (f)|>k.c.  In  £14)  for  N  . (f )  we  set  N  (f)  to  be  any  function  integrating 

o°v^ .  on^A.  if  a  .  =0.  TM  robust  filte^  is 


to 


Rj 


a 


VO  = 


kexp(-iargS^(f)) ,  f  Rj 


(19) 


Ho  »>/<= 


f  eA  -R 


Comments  Note  that  (  17)  and  (  are  the  constraints  on  total  noise  power 
and  on  maximum  deviation  from  the  nominal  signal  for  signals  in  S  ,  respect¬ 
ively.  The  solution  for  S_  and  N  in  any  specific  case  requires  a  numerical 
procedure  for  evaluating  the  constants.  In  principle,  this  procedure  can  be 
an  iterative  procedure  starting  f^om  some  guess  for  and  c.  For  given  k 
and  c  a  solution  for  the  R  and  a^Acan  be  sought  from  (13)  and  (1^).  Then 
the  constraints  (17)  and  CtS  )  are  tested.  A  proof  is  given  in  the  Appendix. 


It  is  also  intersting  to  note  that  as  the  partition  becomes 

finer  the  solution  approaches  that  for  the  "band"  class  fof^N^ln  [4).  In 
the  special  case  where  only  a  noise  power  constraint  is  given  (j=l)  the  ro¬ 
bust  filter  =»  S*/Nj^  =  k^expi-iargS^Cf ) )  on  has  a  constant  magnitude. 

From  (13)  and  (lA),  the  optimum  filter  for  is  not  uniquely  de¬ 

fined  where  is  !<k,c.  On  this  set  the  robust  filter  is  simply 

it  *0*1 

S  (f)/c,  i.e.,  a  filter  matched  to  S  for  white  noise, 
o  o 

IV.  ROBUST  TIME-DELAY  ESTIMATION 


In  a  recent  paper  [8],  the  specified  p-polnt  classes  and  the  band  classes 
of  bounded  spectral  densities  have  been  used  in  obtaining  robust  versions  of 
the  Eckart  filter  [8],  which  is  a  specific  type  of  weighting  function  used  in 
the  computation  of  the  cross-correlation  function  in  time-delay  estimation. 

The  basic  problem  is  to  estimate  the  time  delay  D  between  noisy  signals  re¬ 
ceived  at  two  sensors,  the  sensor  Inputs  being  described  by  v^(t)  s(t) 

+  n.(t)  and  V2(t)  ■  s(t-D)  +  n2(t).  Here  s(t)  is  a  stationary  random  signal 
with  spectrum  S  and  n^^  and  n2  are  uncorrelated  noise  processes  Independent 
of  s,  with  spectra  N.  and  N2,  respectively.  The  technique  of  time-delay  es- 
Clmatlon  considered  In  [8]  is  based  on  forming  an  estimate  of  the  cross- 


spectral  density  G,„(f)  of  the  inputs  and  ^2^  weighting  this  with  a  fi 
ter  function  W,  and'^then  taking  the  inverse  transform  to  obtain  an  estima 


fil- 

imate 


the  cross-correlation  function  R, ^(t).  The  position  of  the  peak  in  R  „(:) 


of  ^ 

is  taken  as  an  estimate  of  the  tinie‘'’^elay  D.  For  low  signal-to-noise  ratio 
at  the  input  one  appropriate  performance  measure  is  the  output  SNR  of  the 

This  is  given  by 

2 


estimate  of 


d(S,Q;W)  = 


_  [.'•W(f)S(f)df]' 


;w^(f)Q(f)df 


(20) 


where  Q(f)  =  Nj^(f)N  (f).  The  spectral  product  Q  and  the  signal  spectrum  can 
both  be  estimated  from  measurements  at  the  system  output.  For  given  S  and 
Q  the  optimum  (Eckart)  filter  maximizing  d(S.Q;W)  is  W  =  S/Q,  and  this  max¬ 
imum  distance  is  d^(S,Q)  =  /S  /Q  df.  ° 

The  results  in  [8]  show  that  for  convex  classes  of  allowable  S  and  Q, 
the  least-favorable  pair  (S^^,  for  (S,Q),  minimizing  d^(S,Q),  and  the 

corresponding  optimum  filter  W  =  form  a  saddle  point  for  d(S,Q;W). 

Thus,  the  saddlepoint  robust  filter  can  be  determined  once  the  least  favor¬ 
able  pair  1^^  known.  Now  since  d  (S,Q)  can  be  written  as  /(S/Q)  Qdf , 

lemma  1  can  be  used  to  obtain  directly  tHe  least  favorable  pair  from  the 
corresponding  hypothesis  testing  problem.  For  S  and  Q  lying  in  power  con¬ 
strained  bounded  p-polnt  classes  the  least  favorable  pair  is  exactly  the 
same  as  for  Wiener  filtering  in  Section  II,  where  Q  plays  the  role  of  a 
noise  spectrum. 

V.  CONCLUSION 

We  have  considered  robust  hypothesis  testing  and  robust  random  signal 
estimation,  deterministic  signal  detection  and  time-delay  estimation,  for 
cases  where  probability  and  spectral  density  functions  belong  to  power- 
constrained  bounded  p-point  classes.  These  classes  are  of  practical  in¬ 
terest  because  the  bounds  can  be  readily  estimated  in  many  situations.  They 
form  a  generalization  of  the  specified  p-point  classes  and  also  of  band- 
classes,  which  can  be  viewed  as  limiting  cast'  of  the  bounded  p-point  classes. 
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APPENDIX: Proof  of  Theorem  in  Section  III. 


Consider  the  denominator  D(N;H)  =  7n|h|  df  in  (12). 

r  2 

D(N;H^)  -  D(Nj^;H^)  =  ^/[k^R  (f)  + - ^ - 

j  J  J  c 


We  have 

R  (f)]  (N-N„)df,  (Al) 
J  K. 


where  R. 
Rj  il 


(f)  =  1  for  f  c  Ri  and  0  otherwise,  and  where  R.  =  A.-R 
0  and  [s  1^/c’  <  k?,  so  J  J  j 

D(N;Hj^)-D(N^;Hj^)  <  ^  k^  /[A^  (f  )N-R^  (f  )N^]df 

<  I  k^[/A  (f)Ndf-v|  ] 
k  >k  J  ^ 


Now  on 


+  I  k^[/A.(f)Ndf-v^.]  +  I  k^[/A.(f)Ndf-Op./k]  (A2) 

k.<k  ^  k.=k  J 

J  J 

The  RHS  of  (A2)  is  0,  so  D(N;H^)  <  D(Nj^;H^)  for  all  N  £  M. 

For  the  numerator,  we  first  show  that  /SH^^df  >  0,  all  S  e  S.  This 
follows  by  noting  that,  as  in  (Al), 


H^(f)  =  +  Rj(f) |S^(f) |/c]exp(-iargS^(f)), 

replacing  S  by  S  +  (S-S  ).  The  integral  containing  S  only  is  independent 
of  S  and  is  lower  boundeE  by  £  /c,  from  (18).  The  second  integral  [integrand 
containing  (S-S  )]  is  real-valued  and,  from  Schwarz's  inequality,  bounded 
by  6  /c  in  abso?ute  value.  Thus  we  consider  minimization  of  this  integral 
which  is 


/[^kjRj(f)  +  R^ (f ) I S^l /c] I S-S^ |cos[arg(S-S^)-argS^]df 

To  minimize  the  cos  term  is  -1,  when  arg  (S-S  )  and  arg  S  have  a  difference 
of  TT.  Without  the  cos  term,  the  integral  is  maximized  for  S  given  by  (13), 
from  a  direct  application  of  Schwarz's  Inequality  and  the  above  observation 
on  phase. 


